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ABSTRACT 

We present measurements of a population of matched radio sources at 1.4 and 5 GHz down to a flux 
limit of 1.5 mJy in 7 deg^ of the NOAO Deep Field South. We find a significant fraction of sources 
with inverted spectral indices that all have 1.4 GHz fluxes less than 10 mJy and are therefore too faint 
to have been detected and included in previous radio source count models that are matched at multiple 
frequencies. Combined with the matched source population at 1.4 and 5 GHz in 1 deg^ in the ATESP 
survey, we update models for the 5 GHz differential number counts and distributions of spectral indices 
in 5 GHz flux bins that can be used to estimate the unresolved point source contribution to the cosmic 
microwave background temperature anisotropics. We find a shallower logarithmic slope in the 5 GHz 
differential counts than in previously published models for fluxes < 100 mJy as well as larger fractions 
of inverted spectral indices at these fluxes. Because the Planck flux limit for resolved sources is larger 
than 100 mJy in all channels, our modified number counts yield at most a 10% change in the predicted 
Poisson contribution to the Planck temperature power spectrum. For a flux cut of 5 mJy with the 
South Pole Telescope and a flux cut of 20 mJy with the Atacama Cosmology Telescope we predict 
a ~30% and ~10% increase, respectively, in the radio source Poisson power in the lowest frequency 
channels of each experiment relative to that predicted by previous models. 

Subject headings: cosmic background radiation - cosmological parameters - Cosmology: observations 
- inflation - surveys 



1. INTRODUCTION 

A main science driver for ongoing high-resolution cos- 
mic microwave background (CMB) experiments is to 
measure the amplitude and tilt of the power spectrum 
of primordial density fluctuations and thereby constrain 
models of cosmic inflation. On large angular scales (mul- 
tipoles £ < 1000), the power spectrum of CMB tempera- 
ture anisotropics is directly related to the primordial den- 
sity power spectrum by well-understood baryon physics. 
However, on smaller angular scales various astrophysical 
foregrounds dominate the measured temperature power 
spectrum and obscure the primordial signal. 

For frequencies less than about 150 GHz a domi- 
nant foreground is the flux from radio-loud galaxies and 
blazars that is unresolved by the CMB experiment beam 
and therefore cannot be masked out (Tegmark & Ef- 
stathiou 1996; Knox 1999; Scott & White 1999; Huf- 
fenberger et al. 2008; Colombo & Pierpaoli 2010; Mil- 
lea et al. 2011). Previous radio surveys have identified 
a significant population of sources from 1 to 30 GHz 
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that have flat or "inverted" spectral energy distributions 
(SEDs) such that their flux is roughly constant or in- 
creasing with increasing frequency (Guerra et al. 2002; 
Colombo & Pierpaoli 2010; Prandoni et al. 2010) and 
are therefore potentially significant foregrounds for CMB 
observations. Physical models for these objects predict 
that they are dominated by synchrotron emission from 
relativistic jets in active galactic nuclei (AGNs; Danese 
et al. 1987; Toffolatti et al. 1998; de Zotti et al. 2005, 
2010; Tucci et al. 2011; Rani et al. 2011) with the large 
spectral indices determined either by dominant emission 
from compact optically thick regions of the jet, early or 
late stages of AGN evolution, or chance variability in the 
source. Such physical models predict that the SEDs of 
these sources should break and decline somewhere be- 
tween tens and hundreds of GHz. The flat and inverted 
radio sources detected around a few GHz should then 
contribute negligible flux at millimeter frequency CMB 
observations (where infrared emission from galactic dust 
becomes a problem instead). 

Thermal synchrotron emission from advection- 
dominated accretion flows (ADAFs) may be also be a 
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significant contribution to the low-frequency radio fore- 
grounds in CMB experiments (Perna & Di Matteo 2000; 
Pierpaoh & Perna 2004; Toffolatti et al. 2005). ADAFs 
may also show inverted spectra up to several tens of GHz 
when the emission region is significantly compact. Using 
the model of Perna & Di Matteo (2000) and considering 
constraints on the source counts from the Ryle telescope, 
ATCA, and Wilkinson Microwave Anisotropy Probe 
(WMAP) surveys, de Zotti et al. (2005) showed that 
the predicted differential number counts of ADAFs at 
30 GHz are at least two orders of magnitude smaller 
than the blazar counts for fluxes greater than ~ 1 mJy 
(see their Figure 14). We therefore only consider SED 
models for blazars with inverted spectra and neglect 
any possible (sub-dominant) contribution from ADAF 
sources. 

Aside from CMB power spectrum measurements, 
radio foregrounds are also a significant foreground 
for Sunyaev-Zel'dovich measurements of galaxy clus- 
ters (Knox et al. 2004; White & Majumdar 2004; de Zotti 
et al. 2005; Lin et al. 2009; Sehgal et al. 2010; Reese 
et al. 2011) and for testing for non-Gaussianity in the 
CMB (Eisner et al. 2010; Curto et al. 2011). 

As CMB experiments obtain better resolution, the 
lower flux limit will decrease for resolving and remov- 
ing radio sources. It is then important to character- 
ize the radio source number counts and spectral indices 
at fluxes below the resolved source flux cut to miti- 
gate the contamination from unresolved sources. While 
WMAP^ was able to remove radio point sources with 
fluxes greater than 0.7 Jy, Planck^ resolves all sources 
with fluxes greater than ~ 0.2 Jy (Planck Collaboration 
et al. 2011). The Atacama Cosmology Telescope (ACT)^ 
and the South Pole Telescope (SPT)'* are able to resolve 
all sources above 20 and 5 mJy respectively (Marriage 
2011; Vieira et al. 2010). Existing models for radio source 
counts and SEDs have largely been calibrated only to a 
flux limit of 30-100 mJy (de Zotti et al. 2005; Tucci et al. 
2011), so the models would have to be extrapolated to 
lower fluxes in order to analyze these ongoing CMB ex- 
periments. 

We present measurements of the population of faint 
radio sources complete to 1.5 mJy at 1.4 and 5 GHz 
collected with the Very Large Array (VLA) in 7 deg^ 
of the NOAO Deep Field South (DFS). Our sources are 
matched at the two observation frequencies allowing us 
to determine the spectral index for each source. To ap- 
ply our source catalog for predicting the unresolved point 
source contamination in the CMB, we compare primar- 
ily with Tucci et al. (2011) (hereafter Til) who present 
models for the number counts and SEDs of radio sources 
calibrated from many different data sets. Our data are 
new primarily in reaching fainter flux densities with spec- 
tral indices of matched sources. 

This paper is organized as follows. We describe our 
observations and the construction of our source catalog 
in Section 2. We derive new fits to the 5 GHz differen- 
tial number counts of our sources and the distributions 
of spectral indices in Section 3.1. We then review SED 
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models that use the fits to the 5 GHz counts and in- 
dices as inputs for extrapolating source fluxes to higher 
frequencies in Section 3.2. We show the impact of our 
new data and models on the predicted high-frequency 
number counts in Section 4.1 and predict the impact on 
measurements of the Planck temperature-temperature 
(TT) power spectrum in Section 4.2 We draw conclu- 
sions about the impact of faint radio sources on future 
CMB measurements in Section 5. In Appendix A we de- 
scribe a statistical model for propagating the errors in 
the measured fluxes and spectral indices into the esti- 
mated high-frequency differential counts and CMB Pois- 
son contribution from unresolved sources. We describe 
the structure of our new source catalog in Appendix B 
and present a sample of the measured source fluxes in 
Table 5. The full catalog is available for download from 
the VizieR database. 

2. SOURCE CATALOG 

We present a catalog of 362 new sources discovered 
with VLA imaging in 7 deg^ of the NOAO DFS at 1.4 
and 5 GHz (see Appendix B for a sample of the catalog). 
The data reduction and catalog construction were per- 
formed with the same methods as in the Faint Images 
of the Radio Sky at Twenty cm (FIRST) survey (White 
et al. 1997). The survey has a resolution of 5 arcsec at 
both frequencies. We include in our catalog only sources 
that are matched in the two bands. There are 700 identi- 
fied sources in the field at 1.4 GHz, so roughly half of the 
sources are matched at 5 GHz to give us our 362 sources. 
There are seven point sources in our 5 GHz catalog that 
are not matched at 1.4 GHz. All seven sources have 
5 GHz fluxes between 1 and 5 mJy, and therefore likely 
have inverted spectral indices such that the 1.4 GHz flux 
is below the FIRST survey detection threshold around 
1 mJy. Lacking spectral index measurements for these 
seven sources, we do not include them in our subsequent 
analysis. However, all our results must be interpreted 
remembering that there are likely more inverted sources 
near our flux limit that would further boost the differen- 
tial number counts of faint point sources. 

The ATESP (Prandoni et al. 2000; Prandoni et al. 
2006) catalog has a similar resolution to ours and has 
118 matched sources in the same frequency bands over 
1 deg^ . The ATESP catalog is complete below 1 mJy and 
is therefore a useful complement to our catalog. We will 
consider both the NOAO DFS and ATESP catalogs in 
tandem when fitting models for the differential number 
counts and SEDs of the 5 GHz source population. 

The observed fluxes and the inferred spectral indices 
(defined by the relation S (X v" where S is the flux, v 
is the frequency, and a is the spectral index) for each of 
our matched sources are shown in Figure 1. There are 28 
and 9 sources with spectral indices ai.4_5 > 0.3 in the 
NOAO DFS and ATESP catalogs respectively, which we 
label as "inverted" spectrum sources (following a com- 
mon convention). All of our inverted spectrum sources 
have 1.4 GHz fluxes < 10 mJy so they would not have 
been identified in earlier surveys with larger flux lim- 
its. The sources with —0.5 < ai 4_5 < 0.3 (129 in the 
NOAO DFS catalog and 55 in the ATESP catalog) are 
labeled as "flat-spectrum" and could also potentially be 
visible at higher frequencies. Again, we see in the top 
panels of Figure 1 that the majority of our fiat-spectrum 



CMB faint inverted radio foregrounds 




-1.50-0.75 0.00 0.75 1.50-1.50-0.75 0.00 0.75 1.50 
a 1.4-5 

Figure 1. Flux in the observation bands of 1.4 and 5 GHz vs. the 
spectral index inferred from these bands. The bottom panels show 
the fluxes extrapolated to 100 GHz using two different SED models 
as described in Section 3.2. The bottom left panel uses a simple 
power-law (PL) model while the bottom right panel uses a model 
from Til that has a "turnover" in the SED at ~ 10 30 GHz. 



sources have 1.4 GHz fluxes below 10 mJy. As a prelim- 
inary illustration of the potential contribution of these 
source populations for high-frequency number counts we 
show the extrapolated fluxes of our sources at 100 GHz 
using two different SED models in the bottom panels 
of Figure 1. We will explain the SED models in the 
next section, but note here that with both a "power- 
law" (PL) SED model and an SED model with a break 
or "turnover" many of our sources are predicted to have 
100 GHz fluxes greater than 1 mJy and will therefore be 
detectable in high-resolution CMB experiments. In our 
fits to the 5 GHz differential number counts in the next 
section, it will also be important to note that there is not 
a strong correlation between the 5 GHz flux and ai.4-5 
in Figure 1. 

In Figure 2 we show the differential number counts at 
1.4 GHz estimated from the NOAO DFS and ATESP cat- 
alogs and compare with the fit to the number counts pub- 
lished by Tegmark & Efstathiou (1996) from the FIRST 
survey at 1.4 GHz (White et al. 1997), which is claimed 
to be complete to 5 = 0.75 mJy Both the NOAO DFS 
and ATESP data match the fit well for S" > 10 mJy but 
the NOAO DFS data clearly give fewer number counts 
for 5 < 3 mJy while the ATESP counts drop below the 
fit for S* < 1 mJy. The incompleteness at low flux in 
both catalogs is primarily caused by the requirement to 
match sources at 5 GHz, which omits steep-spectrum- 
sources at 1.4 GHz that are not detected at 5 GHz. 
This selection effect is apparent in the top left panel 
of Figure 1 where the minimum flux for steep-spectrum 
sources (aiA-5 < —0.5) is much larger than that for 
flat-|-inverted sources. 

3. MODELS FOR SOURCE COUNTS AND 
SPECTRAL ENERGY DISTRIBUTIONS 
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Figure 2. Estimated differential number counts at 1.4 GHz. The 
solid line is the fit to the counts from the VLA FIRST survey at 
1.5 GHz from Tegmark & Efstathiou (1996), which is complete to 
0.75 mJy. 



In this section, we describe the modifications based on 
our data to the Til models for 5 GHz number counts, 
spectral index distributions, and SEDs used for extrap- 
olation of the counts to higher frequencies. The key 
changes to the Til model are from the extension to lower 
5 GHz fluxes. 

3.1. Models for 5 GHz Sources 

We adopt the form of the fit function from Equa- 
tion (1) of Til for the 5 GHz differential number counts 
for flat-|-inverted spectrum sources. 



dN, 



fl-(-inv 



dS-^ 
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Til found best fit parameters no,/ = 47.4 Jy ^sr ^, 
5*0,/ — 1.67 Jy, kf = 0.50 and £f = —0.66. Because of the 
small areas covered by the NOAO DFS and ATESP cat- 
alogs, we do not have good statistics for constraining the 
differential number counts at fluxes greater than about 
100 mJy. We therefore include the differential count es- 
timates (covering 0.1 ^ 5 < 10 Jy) presented in the 
right panel of Figure 2 of Til when deriving new best 
flt parameters for the fltting function in Equation (1). 
In Til, the steep-spectrum number counts are defined 
as the difference between the dN/dS^ model of Toffolatti 
et al. (1998) and Equation (1). Here, we introduce a flt 
function for the steep-spectrum counts that can be used 
to flt our steep-spectrum sources. 
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where no,s = 88 Jy sr 
if = -0.49, m.s = 10.9 Jy-^sr-i, and a^ = 0.33 fit the 
steep-spectrum count model reported in Til. 
When fltting to the differential counts estimated from 
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our catalogs, we multiply Equations (1) and (2) by selec- 
tion functions of the form, 



'/'(S'5;S'min,a) 



1 



-(S5/S„ 



(3) 



to account for incomplete source extraction near the lim- 
iting fluxes in each catalog. We fit separate values of 5'min 
and a for the NOAO DPS and ATESP catalogs. For each 
catalog, we also fit separate S'min and a values for the 
fiat+inverted and steep-spectrum samples because the 
fiux limit of matched sources is a function of the source 
spectral index (see, e.g., the skewed lower bounds of the 
scatter plots in the top panels of Figure 1). 

We fit new values for the flat-|-inverted differential 
number count model in Equation (1) with the following 
procedure. 

1. We create separate histograms for the 5 GHz fluxes 
in the NOAO DFS and ATESP catalogs with 
the optimal histogram bins chosen as prescribed 
by Hogg (2008) (assuming the logarithm of the 
bin widths is constant). The histogram bin val- 
ues with Poisson uncertainties are our estimators 
for the 5 GHz differential number counts. 

2. We combine our number count estimates with those 
presented in the right panel of Figure 2 in Til, 
which compares the 5 GHz count estimates from 
a number of surveys covering the flux range '^ 
100 mJy to ~ 10 Jy. Because our data cover a small 
area, we do not have good estimates of the num- 
ber counts above ^ 100 mJy so the data from Til 
is essential to constrain the "bend" and high-flux 
slope of the flt function in Equation (1). Note that 
we assign new Poisson uncertainties to the data 
points extracted from Figure 2 in Til according 
to the number of sources in each flux bin (derived 
according to the area of each survey). The uncer- 
tainties assigned to our differential count estimates 
and those from Til are therefore consistent. 

3. We extract samples from the joint posterior of all 
the fit parameters in Equation (1) and S'min and a 
(from Equation (3)) for each of the NOAO DFS and 
ATESP data sets using Markov Chain Monte Carlo 
(MCMC) with a likelihood given by the product of 
the Poisson likelihoods for each count estimate in 
each flux bin for each data set (our two catalogs and 
those from Figure 2 in Til). We thin the MCMC 
samples to obtain uncorrelated samples of the joint 
parameter posterior. 

4. We define the best fit parameters as the medians 
of the marginal posteriors for each parameter. 

Using this algorithm, the new flt parameters for Equa- 
tion (1) are, kf = 0.438 ± 0.0142, no = 46.7 ± 
1.49 Jy-^sr-i, £ = -0.755 ± 0.137, and Sq = 2.31 ± 
0.222 Jy. Note in particular that the addition of our 
data favors a much smaller value of kf than the value of 
kf = 0.5 presented in Til, which implies a larger number 
density of faint sources than would be inferred by extrap- 
olating the Til flts to faint fluxes. The logarithm of the 
ratio of likelihoods for our best flt parameters and those 
of Til is 23.6, indicating that our new fit is favored by 
the data with strong significance. The selection function 



minimum flux flts are S'min =1.48 for the NOAO DFS 
catalog and Smin =0.405 for the ATESP catalog. We use 
these flts to determine the minimum flux values when 
creating simulated source catalogs in Section 3.2. 

We adopt a simpler maximum-likelihood procedure to 
fit the steep-spectrum count parameters defined in Equa- 
tion (2) for three reasons: (1) the steep-spectrum sources 
presented in the left panel of Figure 2 in Til have much 
larger uncertainties for constraining the fit function, (2) 
our data turns out to be entirely consistent with the fit 
function parameters derived from the Til results, and 
(3) the steep-spectrum number counts are not important 
for predicting the higher frequency counts at the CMB 
frequencies we are considering in this paper. The bottom 
panel in Figure 3 illustrates how our data are consistent 
with the Til result. 

The distributions of spectral indices in 5 GHz fiux bins 
are shown in Figure 4. We make the simplifying assump- 
tion that the spectral index is statistically independent 
of the 5 GHz flux, which is consistent with Figure 1. 
We again follow Til and fit the histograms with sepa- 
rate truncated distributions for flat-|-inverted (a > —0.5) 
and steep {a < —0.5) spectral types (shown by the red 
solid lines in Figure 4). The dashed black lines in the 
figure show the fits from Til for their lowest-flux bin, 
100 < S/mJy < 158. We find the distributions are 
best fit by a mixture of two truncated Gaussian distri- 
butions, except for the highest-fiux steep-source distri- 
bution, which is well fitted by a single truncated Gaus- 
sian. The parameters of the fits in Figure 4 are given 
in Table 1 along with uncertainties on the flt param- 
eters derived from the 68% conditional uncertainty in- 
terval on the mean of the largest-index Gaussian fit in 
each panel of Figure 4. In detail, we calculate uncertain- 
ties via bootstrap resampling of each catalog, where we 
also sample the number of catalog entries from a Pois- 
son distribution and then sample with replacement the 
catalog entries. We then fit double or single Gaussians 
to each sub-sample shown in Figure 4, rank order the 
fitted Gaussian means, or the larger of the two means 
when two Gaussians are fit, and select the 68% confi- 
dence intervals from the ordered samples. We focus on 
the Gaussian means to understand the possible variation 
(including Poisson uncertainties) in the inverted spectral 
index distributions. 

It is notable for extrapolating fluxes that our best fit 
distributions for the flat+inverted spectrum sources have 
larger tails with positive spectral indices than the Til fits 
in their lowest flux bin. 

Note that while Til consider how source variability 
can change the inferred spectral indices, their modeling 
shows that variability has only a small effect on the spec- 
tral index histograms. We therefore ignore the effect of 
variability here. 

3.2. Flux Extrapolation 

We will primarily use the Til SED models for ex- 
trapolating the 5GHz fluxes to higher frequencies (which 
we will refer to as the "Tucci SED" model). Til con- 
sidered both many different data sets with spectral in- 
dex information as well as physical models for the fre- 
quency dependence of the synchrotron emission from the 
flat-spectrum sources to construct statistical SED mod- 
els that are consistent with observed differential number 
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Figure 3. Observed differential number counts at 5 GHz and fits for the flat+inverted spectrum (top) and steep-spectrum (bottom) 
sources. The red solid lines in each panel are the flts to the NOAO DFS counts (red points and error bars) while the blue dashed lines are 
the fits to the ATESP counts (blue triangles and error bars). The solid and long-dashed black lines show the fits from Tucci et al. (2011) 
for the NOAO DFS and ATESP selection functions, respectively. For all data points, the error bars show the 68% confidence intervals 
assuming the counts in each flux bin are Poisson distributed. 



Table 1 

Fit Parameters for Spectral Index Distributions 



Spectral type 


Flux Bin 
(mJy) 


Mean 1 


Mean 2 


Std. Dev. 1 


Std. Dev. 2 


Weight 


flat-l-inverted 
flat+inverted 

steep 

steep 
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Figure 4. Distribution of spectral indices for flat+inverted (left) 
and steep-spectrum (right) sources in two 5 GHz flux bins (top 
and bottom rows) measured from the combined catalogs of the 
NOAO DFS and ATESP fields. The flux bin ranges in the row 
side panels are in mJy. The black dashed lines show the flts from 
Til in their [100, 158] mJy flux bin, which is the lowest flux bin 
they consider. 



counts at frequencies from 5 to several hundred GHz. 
We will show in this section that our modified fits to the 
5GHz flat+inverted differential counts and the spectral 
index distributions for fluxes less than 100 mJy further 
improve the fits of the extrapolated differential number 
counts at low fluxes with external data sets. 

The Tucci SED model assumes that steep-spectrum 
sources follow a PL SED {S (x v"^) with spectral indices 
ai.4_5 — Aa and Aa Gaussian distributed to model ob- 
served steepening at higher frequencies (see Equation (9) 
in Til). In the description of the steep-spectrum SED, 
Til state that, "A small percentage of flattening or 
upturning sources is also included." We find that we 
can reproduce their plots of extrapolated steep-spectrum 
source counts by assigning 2% of our steep-spectrum 
sources to have spectral indices drawn from a Gaussian 
distribution with mean —0.3 and standard deviation 0.2. 

We adopt the "C2Ex" model from Til for flat- 
spectrum sources (—0.5 < ai.4_5 < 0.3) that deter- 
mines the distribution of break frequencies according to 
a physical model of the size of the optically thick syn- 
chrotron emitting region in FSRQ and BL Lac sources. 
The model requires knowledge of the redshift distribu- 
tions of the sources of each type as well as the dis- 
tributions of Doppler factors that adjust the predicted 
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flux from a fiomogeneous spherically symmetric model 
to model an asymmetric jet. Our sample of faint flat- 
spectrum sources almost certainly has a different red- 
shift distribution than that assumed by Til for their 
high-frequency predictions. This is because for a homo- 
geneous source sample, fainter fluxes imply higher red- 
shifts, while at fixed redshift fainter fiuxes imply sources 
at different evolutionary stage that will again have dif- 
ferent redshift distributions. However, in the absence of 
further information about our catalogs, we simply use the 
same distributions for the flat-spectrum break frequen- 
cies as Til (shown in Figure 7 of Til). This is probably 
not a terrible approximation because the distributions of 
the break frequencies span several orders of magnitude 
with the mean following a simple scaling with 5 GHz flux 
and the width determined largely by the assumed size of 
the emitting region. So modifications in the assumed 
redshift distribution would yield only minor corrections 
in the final predicted number counts. 

For inverted spectrum sources {aiA-5 > 0.3) Til use 
the model. 



s{,y,e) = So 



where 



I ~ —A {I — Iq) exp 



1-e" 



2 



hY 



(4) 



(5) 



with k = -0.1 and ui = 0.53 (see Figure 9 in Til). Note 
that / is therefore defined with support [— oo, /q) and ki 
must be less than 0. The parameter ui^ in Equation (4) 
is determined from the distribution of peak frequencies 
in Tabic 3 of Til. Again we adopt the distribution for I 
from Til without modification. 

The Til SED model as applied in this paper has two 
main drawbacks. First, the model requires specification 
of both the distributions of "break" frequencies where 
the SEDs of inverted and flat-spectrum sources turn over 
as the radio emission probes different physical regions of 
the source and of the distributions of the spectral indices 
after the break. Second, and somewhat related, the Til 
SED model has a large number of parameters with un- 
specified uncertainties so that it is difficult to quantify 
our confidence in the extrapolated number counts pre- 
dicted with this model. For comparison then, we also 
consider the physically unrealistic but simple PL SED 
model, 



S{u;Sfi,a) ^ So 



^0 



(6) 



with vo — ?) GHz. The PL model will overpredict the 
number counts at frequencies larger than the typical 
break frequencies for a given source type. But, the PL 
model has the advantage that we can easily propagate 
our uncertainty in the SED model parameters, Sq and 
a, into the extrapolated number count predictions and 
therefore gain some understanding about the relative un- 
certainties in the number counts from the errors in the 
data versus the SED model. 

We propagate the uncertainty in the PL SED model 
by computing the marginal posterior probability distri- 
bution for the extrapolated flux of each source marginal- 



izing over the SED parameters, 5*0 and a, and incorpo- 
rating observational uncertainties in the flux likelihood 
model. We describe the details of this extrapolation 
method in Appendix A. The marginal posterior prob- 
ability for the extrapolated flux allows for full uncer- 
tainty propagation and therefore robust foreground re- 
moval (conditioned on the choice of SED model). Such a 
procedure becomes less useful when the available source 
catalog is an incomplete sample of the population of 
sources contributing to the foregrounds and when there 
are many SED parameters to marginalize over leading to 
large marginal uncertainties. Our catalogs unfortunately 
meet the former condition because of the relatively small 
fraction of the sky covered (leading to a dearth of high- 
flux sources observed) and the Tucci SED model likely 
meets the latter condition. 

Putting aside the incompleteness of our catalogs for 
fluxes > 100 mjy, we predict the high-frequency number 
counts with the PL SED model using the complete error 
propagation methods of Appendix A. Assuming a log- 
Normal likelihood for the observed fluxes, the marginal 
posteriors for the extrapolated fluxes can be computed 
analytically (and therefore quickly). The resulting ex- 
trapolated flux confidence intervals can be used as indi- 
cators of the extrapolation error when the SED model is 
well-constrained. 

However, with the many parameters in the Tucci SED 
model that, upon marginalization, further contribute to 
the extrapolated fiux uncertainties, we instead use the 
Til flux extrapolation method to study the impact of 
our new data on the "central" extrapolated count pre- 
dictions. In Til the differential number counts at high 
frequencies are predicted by first simulating a source cat- 
alog by sampling sources from the fits to the 5 GHz differ- 
ential counts and sampling spectral indices from the fits 
to the 5 GHz spectral index histograms, and then plug- 
ging the sampled values into the SED models for each 
spectral class. The extrapolated simulated fluxes can 
be histogrammed to estimate the high-frequency number 
counts in the same way that the observed low-frequency 
number counts were estimated. Because the predictions 
from any one simulated catalog can be subject to random 
fluctuations (mostly at high fluxes where the counts are 
low) we generate predictions for the Tucci SED model 
by simulating 10 catalogs with 10^ sources each and then 
computing the medians of any summary statistics of the 
extrapolated catalogs. 

4. HIGH-FREQUENCY PREDICTIONS 

We now apply the flux extrapolation methods de- 
scribed in Section 3.2 to predict the high-frequency dif- 
ferential number counts and Poisson contributions to the 
CMB power spectrum. We consider two sets of simulated 
catalogs for generating predictions with the Tucci SED 
model. First we simulate catalogs with exactly the same 
parameters as in Til (labeled "Tucci et al. 2011" in the 
flgures) but with a lower fiux limit of 0.4 mJy, which is 
where our fit to the ATESP selection function (see Equa- 
tion (3)) has a value of 0.5. To reach this flux limit we un- 
fairly apply the spectral index distributions from Til for 
the fiux bin [100,158] mJy to all lower fiuxes. Our second 
set of catalogs (labeled "This work" in the figures) uses 
the modified dN/dS fit with kf = 0.438 for fiat -(-inverted 
sources as found in Figure 3 and the spectral index distri- 
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butions fit in Figure 4 to simulate the spectral indices of 
sources with fluxes lower than the lower limit of 100 mJy 
in Til. As described in Section 3.2, we use the Til 
models for both our simulated catalogs for the flat and 
inverted spectrum break frequencies and high-frequency 
spectral indices. 

4.1. Extrapolated Number Counts 

To validate the differential counts model based on our 
catalogs we first compare with the measured differential 
counts at 15.7 GHz in the 9C and IOC surveys (AMI 
Consortium et al. 2011) in Figure 5. The combined IOC 
survey covers 12 deg'^ to a flux limit of 0.5 mJy while the 
9C survey data cover several hundred square degrees to 
a limit of 10 mJy. While this provides excellent statistics 
for the estimates of the differential number counts, the 
surveys only have a resolution of 30 arcsec compared to 
our catalogs that have a resolution of 5 arcsec. It is there- 
fore possible that the number of sources observed with 
30 arcsec resolution is artificially reduced relative to our 
model because resolved sources at 5 arcsec are merged to- 
gether when they happen to fall close together in projec- 
tion on the sky. Out of the 362 sources in our NOAO DFS 
catalog, we find 29 source pairs that lie within 30 arcsec 
of each other and 10 source triplets within a circle of 
30 arcsec diameter. To predict the 15.7 GHz differen- 
tial counts for 30 arcsec resolution we randomly choose 
sources to merge in our simulated catalogs so that the 
fractions of merged pairs and triplets match those in the 
NOAO DFS catalog. Because of low statistics we have 
made no attempt to choose sources to merge based on 
any other criteria than random selection. The resulting 
predicted number counts are shown by the solid blue line 
labeled "This work, low-res" in Figure 5. Our model is 
a much better match to the IOC survey number counts 
than either the Til model or the de Zotti et al. (2005) 
model for fluxes less than 8 mJy. The better fit at faint 
fluxes is mostly due to the increased 5 GHz counts from 
the fit we found in Figure 3. For the fluxes from 10 mJy 
to 60 mJy the 9C survey number counts fall well below all 
of the model predictions, while all the models are in good 
agreement with the observed counts for larger fluxes. 

In Figure 6 we recreate Figure 12 from Til to compare 
our predicted differential number counts at 30 GHz with 
several other surveys. The steep-spectrum sources are 
a small but non-negligible contribution to the number 
counts at 30 GHz. The new fit to the low-flux slope of 
the 5 GHz differential counts in Figure 3 and the larger 
values of the spectral indices shown in Figure 4 lead to a 
significant difference in the total 30 GHz counts for fluxes 
between 1 and '^ 50 mJy. Our fit leads to greater con- 
sistency than the Til model with the counts measured 
by the Cosmic Background Imager (CBI) (Mason et al. 
2003) and the Sunyaev-Zel'dovich Array (SZA) (Mu- 
chovej et al. 2010), but is in greater disagreement with 
the Green Bank Telescope (GBT) survey by Mason et al. 
(2009) (all of which were previously presented and com- 
pared in Til). The disagreement between our predic- 
tions and the GBT survey measurements might be ex- 
plained in part by the lower resolution of 24 arcsec for 
the GBT survey as well as the targeted source selection 
from the NVSS. For all fiuxes less than ~ 0.1 Jy our 
predicted counts are now larger than the de Zotti et al. 
(2005) model (shown by the dotted green fine). Also 



^ 10' - 









^^ 










^(< 






1 










J 


-r 










^j^ 











\ J* 


'^M 






^/' 


', ~ 








^/ 






























/ 













Model 






Tuccietal. (2011) 




This work, low-res 






TUs work 





De Zotti 


Area [sq. deg.] 

• 11.96 
A 26.86 


■ 


47.23 


+ 


47.83 


S 


124.6 


* 


520 



10"' 10- 10"' 

S(v= 15.7GHz) [Jy] 



10" 



Figure 5. Comparison of our predicted counts at 15.7 GHz with 
those measured in the 9C and IOC surveys by AMI Consortium 
et al. (2011). The dash-dotted green hne shows our model pre- 
diction while the solid blue line shows the same model prediction 
after merging sources to be consistent with the 30 arcsecond reso- 
lution of the 9C and IOC surveys. In both cases we draw sources 
from our fits to the 5 GHz counts and 1.4-5 GHz spectral index 
distributions and then extrapolate the fluxes to 15.7 GHz (before 
merging the sources for the solid blue line). The dashed red line 
instead uses the 5 GHz counts and index distribution fits provided 
in Til (unfairly applied at fluxes lower than those considered in 
Til). The long-dashed purple line shows the model from de Zotti 
et al. (2005). The data points are from Table 6 of AMI Consortium 
et al. (2011) with Poisson error bars added. 
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S(v = 30GHz) [Jy] 

Figure 6. Predicted differential number counts at 30 GHz using 
the Til SED model to extrapolate the modeled 5 GHz counts. 
The thin lines show the steep source counts (at the bottom of the 
figure), the slightly thicker lines show the flat-|-invertcd counts, 
and the thickest lines show the sum of the steep and flat -f inverted 
counts. The solid blue lines use our fits to the 5 GHz counts and 
1.4—5 GHz spectral index distributions while the dashed red lines 
use those from Til (as in Figure ??). The dotted green line is 
the model from de Zotti et al. (2005). The shaded bands showing 
measurements at 33 GHz are copied from Figure 12 in Til. 



note that our predicted number counts are significantly 
larger than those of Colombo & Pierpaoli (2010) (their 
Figure 5) from a simulation based on the NVSS source 
catalog. 
Finally, ACT (Marriage 2011) and SPT (Vieira et al. 
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Predicted Poisson Power for Planek Lower-frequency Channels 
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Figure 7. Predicted differential number counts at 148 GHz with 
the same models as in Figures ?? and 6. The dark gray shaded band 
denotes the range of extrapolated number count predictions when 
the faint slope, kf, for the fiat+inverted 5 GHz number counts is 
varied over its 68% confidence interval. The data points are taken 
from Marriage (2011) for ACT and Vieira et al. (2010) for SPT. 
The light gray shaded band shows the 95% confidence intervals 
for the number counts when each source flux in the NOAO DFS 
and ATESP catalogs is extrapolated using a power-law SED as de- 
scribed in Appendix A. The power-law SED band is intended to 
show the range of extrapolated count uncertainties when marginal- 
izing over SED parameters given only the NOAO DFS and ATESP 
catalogs. 



2010) have recently presented catalogs of resolved point 
sources at 148 GHz, which we compare with our model 
in Figure 7. Again, our model predicts larger number 
counts than the Til and De Zotti models for fluxes less 
than ~ 50 mJy. The uncertainty in the predicted counts 
from the 68% uncertainty range in fc/ is shown by the 
gray shaded band in Figure 7 and is much less than the 
difference between models for S < 50 mJy. All the mod- 
els are consistent with the SPT measurements but the 
increased counts for faint fluxes predicted by our model 
is in tension with the lowest-flux ACT measurement. The 
model predictions at 148 GHz depend sensitively on the 
SED model parameters, in particular the flat-spectrum 
source break frequency distribution described in Til and 
Section 3.2. So it is possible that the tension between 
our model prediction and the ACT measurement indi- 
cates that the break frequencies are really smaller than 
we have modeled for the flat-spectrum sources in our 
catalog with 5 GHz fluxes less than 10 mJy. However, 
the break frequencies are degenerate with the spectral 
indices after the break in the SED when comparing with 
high-frequency number counts so we cannot make un- 
ambiguous conclusions about the SED model from the 
comparison in Figure 7. 

We have also plotted the 95% confidence intervals for 
the extrapolated number counts as a gray shaded band 
in Figure 7 assuming a PL SED model and using the er- 
ror propagation methods described in Appendix A. The 
confidence intervals for the PL SED model include both 
uncertainties in the SED parameters and the Poisson un- 
certainties from the limited size of the NOAO DFS and 
ATESP catalogs (the PL SED predictions use only the 
NOAO DFS and ATESP catalogs as input). When ex- 



Preq. 
(GHz) 



'-'cut 

(Jy) 



C'i Tucci 



Ct This Work 
(mK2) 



— — Tucci et al. 


(2011) 


30 


0.48 


2.8E-02 


2.7E-02 


This work 




44 


0.58 


6.6E-03 


6.2E-03 


■ - ■ De Zotti 




70 


0.48 


8.2E-04 


7.7E-04 






100 


0.34 


1.6E-04 


1.5E-04 


Data set 




143 


0.21 


3.3E-05 


3.2E-05 


-•- ACT 




217 


0.18 


1.5E-05 


1.5E-05 


-*-,SPT 

















trapolating source fluxes over such a broad frequency 
range the PL SED is of course a poor model, but we 
include it in Figure 7 to demonstrate the expected mod- 
eling uncertainties for the flux extrapolation. A simi- 
lar error propagation for all the SED parameters in the 
Tucci model could yield even larger uncertainties. But 
the Tucci SED model is calibrated against many other 
datasets so the priors on the SED parameters might limit 
the increase in the uncertainties relative to the PL SED 
model when all SED parameters are marginalized. We 
have left this investigation for future work. 

4.2. CMB Power Spectrum 

As shown in, e.g., Tegmark & Efstathiou (1996); Scott 
& White (1999) the Poisson contribution to the CMB 
power spectrum from unresolved point sources is 



C,{v)^g\v) 



Sou 






(7) 



where g{v) converts from intensity at frequency v to 
temperature in /iK, and S'cut is the minimum flux (in 
a given channel) at which point sources can be resolved 
and masked out or otherwise removed. The same radio 
sources also contribute a clustering foreground to the 
CMB power spectrum, but this is expected to be sub- 
dominant to the Poisson power for the source intensities 
we are considering here (Scott & White 1999). 

In Table 2 we evaluate Equation (7) for the lower fre- 
quency Planck channels using the 5cut values from the 
Early Release Compact Source catalog (ERCSC) given 
in Table 3 of Planck Collaboration et al. (2011). We 
have omitted predictions for the Poisson power for the 
three highest frequency Planck channels both because 
the unresolved radio sources become a sub-dominant 
foreground at those frequencies and because the flux ex- 
trapolation becomes increasingly unreliable. 

We show similar predictions for the Poisson power for 
ACT and SPT in Table 3 assuming S'cut = 0.005 Jy 
for SPT (Vieira et al. 2010) and S'cut = 0.02 Jy for 
ACT (Marriage 2011). 

We further compare the different model predictions for 
the Poisson power in Planck, ACT and SPT channels in 
Figure 8. Each panel shows the predicted Poisson power 
for our model and the De Zotti model, both normalized to 
the prediction from the Til model as a functions of Scut- 
The top six panels show the predictions for the six lowest- 
frequency Planck channels while the bottom two panels 
show predictions for the lower ACT and SPT frequencies 
and Scut ranges. We have plotted only those frequency 
channels that have predicted radio source Poisson power 
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Table 3 

Predicted Poisson Power for ACT and SPT 



Table 4 

Scalar Spectral Index Bias from Wrong Point-source Model with Four 
Low-frequency Planck Channels. 



Freq. 
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Ci This Work 
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0.005 
0.005 
0.02 
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0.02 
0.02 



2.5E-06 
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2.7E-06 
3.2E-07 
1.5E-06 
1.6E-06 



3.1E-06 
7.6E-07 
3.1E-06 
4.0E-07 
1.7E-06 
1.8E-06 



much larger or comparable to the infrared source Pois- 
son power as predicted by the model from Equation (4) 
of Millea et al. (2011) (using their amplitude normal- 
ized to Planck and SPT). We plot the predicted infrared 
source Poisson power with black crosses or arrows (when 
outside the scale of the panel) only where the predic- 
tion is comparable to or larger than the radio source 
Poisson power (i.e. the highest frequency panel for each 
S'cut range). The circular points and red lines in Figure 8 
show the predictions using the mean of the marginal pos- 
terior for kf from Equation (1) while the blue triangles 
and lines show the model from de Zotti et al. (2005). 
The red shaded bands show the range of predictions for 
kf values spanning the 68% confidence intervals of the 
marginal posterior given the 5 GHz number count ob- 
servations. The gray shaded bands are similar, but also 
include the variation of the 1.4-5 GHz spectral indices 
of the flat+inverted sources over their 68% confidence 
intervals as reported in Table 1. For Planck, our model 
predicts Poisson power that matches the prediction from 
the Til model to within 10% for all S'cut values. This is 
because our data modify the Til model at 5 GHz only for 
source fluxes < 100 mJy while the Planck Poisson power 
is mostly determined by sources with fluxes just below 
S'cut > 100 mJy. The De Zotti model predicts higher 
Poisson power than the Til model and our model be- 
cause of an excess of sources with fluxes just below S'cut- 
This is the reason for the different dependence on Scut of 
the two models in the 70 and 100 GHz panels of Figure 8. 
At the ACT and SPT frequencies in the bottom 2 panels 
of Figure 8, our model predictions are consistently above 
both the De Zotti and Til models, as also indicated in 
Figure 7. The excess faint sources in our catalogs predict 
increased Poisson power in ACT and SPT of 5%-30% rel- 
ative to the Til model, ignoring any uncertainties in the 
Tucci SED model parameters. 

In Figure 9, we plot the different Poisson power model 
predictions in the four lowest-frequency Planck channels 
compared with the TT CMB power spectrum (predicted 
with Emu CMB Schneider et al. 2011) and the expected 
errors from the average noise per pixel and beam smear- 
ing. Figure 9 makes it clear that the Poisson contri- 
bution to the total observed power is comparable to or 
sub-dominant to the noise in each channel. While the 
resolution is better in higher- frequency Planck channels, 
the Poisson power from millimeter sources will dominate 
over the radio sources we are modeling. We also plot in 
Figure 9 the 95% confidence intervals on the predicted 
Poisson power assuming a PL SED derived using the 
methods in Appendix A. The sizes of the uncertainties 
for the PL SED are comparable to the variation in the 
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Poisson power for the range of Scut values in the ERCSC 
and arc much larger than the variation between the Tucci 
SED models for a fixed Scut value. 

We showed in the upper panels of Figure 8 that our 
5 GHz number count measurements predict no more 
than 10% changes in the predicted radio source Poisson 
power in any Planck channel relative to earlier models. 
We now consider whether our data yield 5 GHz count 
measurements that provide a sufficient model to remove 
the bias in constraints on the scalar spectral index, n^, 
from radio foregrounds in the Planck temperature power 
spectrum. We adopt the Fisher matrix formalism as de- 
scribed in Huterer & Takada (2005) to predict the bias 
on Ug assuming our model with the 5 GHz number count 
slope kf ~ 0.452 (our +la fit value) is the "truth," but 
that our model with kf = 0.423 (our —la fit value) is 
used to subtract the Poisson power when analyzing the 
data. For computing the Fisher matrix we use the power 
spectrum noise model including the Poisson contribution 
from White (1998) as well as the optimal linear estima- 
tor therein for the CMB power spectrum given the five 
lowest-frequency Planck channels. The inferred biases 
on Us assuming the same Scut value for each channel are 
shown in Table 4. The last column of the table also shows 
the biases normalized by the marginal Fisher matrix 1-cr 
uncertainties on Hs- We marginalized over 

five other cosmological parameters (as in Schneider et al. 
2011). For all Scut values, the biases on Ug are smaller 
than the la marginal Fisher errors. So, by this measure, 
the remaining uncertainties in the 5 GHz counts fit from 
our catalogs are not significant for Planck parameter es- 
timation (neglecting uncertainties in the SED models). 

5. DISCUSSION 

We have shown that a previously unconsidered popula- 
tion of faint radio sources at 1.4 and 5 GHz with flat and 
inverted spectral indices could significantly contribute to 
the unresolved point source contamination in the power 
spectrum measurements in ongoing high resolution CMB 
experiments. However, the measured differential num- 
ber counts of faint radio sources are not large enough to 
significantly modify previous predictions for the Poisson 
power in the Planck temperature power spectrum. Our 
source catalog is unique in the combination of flux sen- 
sitivity to 1.5 mJy, resolution of 5 arcsec, area covered 
of 7 deg^ and matching of sources in the 1.4 and 5 GHz 
frequency bands. 

We found that existing fits to the differential num- 
ber counts of 5 GHz flat and inverted spectrum sources 
must have a shallower slope for fluxes less ~ 100 mJy 
in order to be consistent with our catalog. We also 
found that the distribution of 1.4-5 GHz spectral in- 
dices of sources with 5 GHz fluxes less than 100 mJy 
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Figure 8. Predicted CMB Poisson power from unresolved radio point sources as a function of the minimum flux, Scut, at wliicli point 
sources can be resolved and excised. The model from this paper ("This work") and the model from de Zotti et al. (2005) are shown by the 
red circles and blue triangles, respectively, normalized by the model from Til. The red shaded bands with the dashed outlines show the 
range of the predicted Poisson power when the flat+inverted 5 GHz differential number count slope kj (see Equation (1)) is varied over its 
68% confidence interval (when fit to our data and that in Til). The gray shaded bands show the range of predicted Poisson power including 
both the aforementioned variation in kf as well as the variation of the 1.4-5 GHz spectral indices over the 68% conditional confidence 
intervals in the mean of the inverted spectral index distributions (see Figure 4 and Table 1). At frequencies where the Poisson power from 
infrared sources is comparable or larger than that from radio sources we have included either a black cross or black arrows (when the Til 
normalized values do not fit within the plot range). In both the top and bottom panels, the millimeter source Poisson power dominates 
the radio source power for higher frequencies than those shown but is highly sub-dominant for all but the highest shown frequency panels. 



is skewed toward larger values than the distributions for 
higher flux sources. Taken together these changes to the 
5 GHz source counts imply increased differential counts 
when the fluxes are extrapolated to higher frequencies 
and therefore increases in the predicted level of Poisson 
power from unresolved point sources in CMB measure- 
ments. 

The quantitative predictions of high-frequency number 
counts depend on the choice of SED model for extrapolat- 
ing the measured 5 GHz fluxes. And the uncertainties in 



the predicted high-frequency counts depend on the prior 
constraints on SED model parameters given a choice of 
SED model. Our main results assume a physically mo- 
tivated SED model (from Tucci et al. 2011) but we also 
investigated a simple PL SED model to assess the depen- 
dence of our predictions on SED modeling uncertainties. 
In general we found that propagating the uncertainties 
from both the limited size of our catalog and two un- 
known SED model parameters can lead to high frequency 
prediction uncertainties that are much larger than the 
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Figure 9. CMB power spectra in tlie four lowest-frequency Planck channels. The Poisson contribution to the power from unresolved point 
sources is shown for three different values of Scut using both the updated 5 GHz number count model and spectral index distributions 
("This work") and the model from de Zotti et al. (2005). The line widths for our model ("This work") denote the prediction range when kf 
is varied over its 68% confidence interval. The Tf f model predictions are always within the line widths of our model. The black line and 
gray shading show the TT CMB power spectrum and expected errors from pixel noise and beam smearing. The orange shaded band shows 
the 95% credible intervals for the Poisson contribution from the sources in the NOAO DPS catalog and assuming a power-law ("PL") SED 
model with Scut = 0.4 Jy. 



differences between the predictions of competing SED 
models. However, we have argued that systematic offsets 
in the high-frequency predictions may persist even when 
the prediction uncertainties are large so that the pre- 
dictions calibrated against our catalog will be important 
considerations for analyzing CMB experiments. As one 
example we find that if the four lowest-frequency Planck 
channels are optimally combined to constrain Ug from 
the CMB temperature power spectrum, the constraints 
are biased by less than the size of the la uncertainties 
if the "wrong" point source number count model is used 
to subtract the Poisson power. We limited our forecasts 
in this example to the five lowest frequency Planck chan- 
nels to avoid modeling other more dominant foregrounds 
at higher frequencies. So, our forecasts are simply an 
illustration of the significance of the difference between 
point source models rather than predictions of the opti- 
mal Planck constraints when all channels are considered, 
but we argue that our catalog shows that radio source 
models at 5 GHz are sufficient for removing Planck Pois- 
son power as a systematic. 

We outlined methods in Appendix A for propagating 
both measurement and SED model uncertainties when 
extrapolating fluxes on a source-by-source basis. A rig- 
orous approach to predicting the Poisson power from un- 
resolved point sources in CMB observations could apply 
the methods of Appendix A to a combined catalog of all 
available radio survey data at relevant frequencies. We 
expect that significant constraints on the uncertain SED 
model parameters could be imposed by such a joint anal- 
ysis of available data, thereby removing a key uncertainty 
from the present work while also providing useful physi- 



cal parameterizations of the radio Poisson foreground (as 
used in Paoletti et al. 2011). However, we leave this for 
future investigations. 
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APPENDIX 

ERROR PROPAGATION FOR 
EXTRAPOLATED FLUXES 

For a given SED model S{iy, 6) with parameters 0, we 
assume a log-normal likelihood for the observed fluxes of 
the zth object in the catalog. We choose a log-normal 
likelihood because it has strictly positive support and 
becomes nearly symmetric about the mean when the flux 
errors are small or the flux is large. 

To propagate our uncertainty in the measured fluxes 
and SED parameters to the extrapolated fluxes, we com- 
pute the marginal posterior distribution for the extrapo- 
lated flux, 



P{Su\Si^ 



dep{s,\e)p{e\s^). 
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The first term gives a deha function, P{S^\0) = 5d{S„ — 
S{i',0)), which fixes one of the SED parameters as a 
function of the remaining SED parameters and the given 
value of S^. For a PL SED, the delta function effectively 

sets a = ln(|^)/ln(^) (or 5*0 = So{S^,a)). The 

marginalization in Equation (Al) then becomes 



PiS,\S,)o, 



1 



S^ln{iy/iyQ) 



dSoPiS,\SQ,aiSo,S,))P{So), 

(A2) 

where the coefficient in front of the integral comes from 
the transformation of variables in the delta function. We 
assume a conjugate log-normal prior for 5*0 with mean 
^iSo s-nd standard deviation asg. 

Performing the integration over Sq, the final expression 
for the marginal posterior for the extrapolated fiux of an 
individual source for a PL SED is 
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where, if we assume vq — 1.4 GHz, 
_ ln(5/z/o) 



X5 



\n{v/vQ) ' 
and with cti^s = as^^^^J 81,4.^5, 
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and, 

As^ = j2TTcr'g^ exp 



-(1 + asj 2/^s„ + (1 + as^)al^ 



We then estimate the mean differential number counts 
by summing over the (normalized) posteriors in bins in 
flux, 
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where ASk = S^axM - 5'inin,fc and fc = 1, . . . , TVbins index 
bins in 5*. 
The Poisson contribution to the CMB power spectrum 



Table 5 

Random Sub-sample of the Catalogue of Matched Sourecs in the NOAO 

DFS. 
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30.86276 


-4.98161 


3.90 


0.27 


7.30 


0.15 


-0.49 


33.28917 


-4.09886 


0.98 


0.35 


6.11 


0.14 


-1.44 


33.34937 


-4.33320 


13.54 


0.38 


65.49 


0.17 


-1.24 


33.88087 


-4.68231 


3.92 


0.46 


7.18 


0.15 


-0.48 


31.76899 


-3.81961 


10.64 


0.30 


34.55 


0.15 


-0.93 


33.44730 


-3.71314 


10.47 


0.30 


33.00 


0.14 


-0.90 


33.56511 


-5.52509 


3.39 


0.58 


4.84 


0.15 


-0.28 


31.06313 


-4.72039 


11.24 


0.30 


30.47 


0.17 


-0.78 


33.19339 


-3.61151 


1.97 


0.49 


18.54 


0.15 


-1.76 


32.27030 


-5.10399 


1.07 


0.28 


4.04 


0.14 


-1.04 


31.32263 


-5.27310 


1.69 


0.38 


2.24 


0.16 


-0.22 


32.76205 


-4.89393 


4.67 


0.36 


2.12 


0.15 


0.62 


34.04510 


-4.43530 


21.66 


0.28 


23.80 


0.15 


-0.07 



can be written as, 
i=l -^0 

JVso„rcc= 

- E 



dSS^PiS,\S,) 



1=1 



xw-^exp 



X 1 - Erf 



(as, -f 3) ((as„+ 3)4^+2^*5 J 



(as„ -t- 3)cr|_^ 4- ^is,, - ln(S'c 



(AlO) 



where the index i in the final equality is implicit in all 
the flux posterior parameters and "Erf" denotes the error 
function. 

NOAO DFS CATALGOUE 

A random sub-sample of the catalog of 362 sources in 
the NOAO DFS matched at 1.4 and 5 GHz is shown in 
Table 5. The spectral index values in the flnal column are 
derived from Columns 3 and 5 via the relation aiA-5 = 
log(S'5/5i.4)/log(5/1.4). The full catalog is available for 
download from the VizieR database. 
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